Hubbard model with SU (4) symmetry 
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In the model considered, the nonlocal interaction of the fermions in different sublattices of a 
bipartite lattice is introduced. It can also be regarded as local interaction of fermions with opposite 
"hypercharge" . The corresponding term in the Hamiltonian is SU (4)-invariant and appears to 
be the most tractable version of the SO (5)-invariant model that unifies antiferromagnetic and 
d-wave superconducting order parameters. The model has been studied primarily in the weak 
interaction limit and in the mean field approximation. Near the half-filling the antiferromagnetic 
critical temperature has a peak. However, the superconducting transition takes place when the Fermi 
surface crosses the area where the density of states is of order of inverse coupling coefficient. Thus, 
in mean-field approximation, there exist an interval of values of the chemical potential, for which the 
system is a superconductor for arbitrary high temperatures. The temperature dependence of specific 
heat, Hall coefficient, and DC conductivity in the normal phase agrees with that experimentally 
observed in high-T c cuprates. 

PACS numbers: 71.10.Fd, 74.20.Mn, 74.25.-q 



I. INTRODUCTION 

Among the properties of the high-T c cuprates that have 
long defied explanation are the proximity of the anti- 
ferromagnetic (AF) and d-wave superconducting (dSC) 
phases below the critical temperature and the anomalous 
temperature dependence of kinetic and thermodynamic 
quantities above the transition point. A recently pro- 
posed concept of an SO (5) symmetry between AF and 
dSC phasestl aims to explain the former as well as the 
resonance mode observed in. spin- flip neutron scattering 
on YBCO.I3 Several groupsHTj have constructed micro- 
scopic models with exact SO (5) symmetry, and it has 
been arguedQ that the two-dimensional (2D) Hubbard 
model has approximate SO (5) symmetry. 

The properties of the microscopic models with sym- 
metries higher than SO (5) were also investigated earlier. 
In Refs. [?]]| the large-n limit of SU (n) model has been 
studied and 1 fn expansion has been applied. It has been 
found that in the strong coupling limit the ground state 
breaks translational symmetry and represents a density 
wave in which each site forms a dimer with one of its 
nearest neighbors. As the doping increases, a "kite" state 
with charge-density wave and no charge gaps forms. In 
the weak interaction limit the "flux" state with full trans- 
lational symmetry and gap vanishing at discrete points 
in momentum space was predicted. However, it was 
shown that at large n the ground state does not have 
off-diagonal long range order. 

In Ref. [)]an SO (6) model has been suggested, in which 
AF, dSC, and flux phase are unified. This and the subse- 
quent workEJ have shown that the pinning of the Fermi 
level near a Van Rave singularity can explain the ob- 
served stripe phasesO in cuprate superconductors. 

This paper presents a study of a relatively simple two- 
leg ladder model that possesses SU (4) symmetry, which 
is higher than SO (5) symmetry. The introduced inter- 



action favors the state, in which each rung is occupied 
by the pairs of fermions corresponding to the opposite 
legs of the ladder, or "flavors" . This model is associated 
directly to a 2D bipartite lattice, since each rung can be 
related to a site in one of the sublattices. Then one leg 
will be simply this sublattice, while the other leg will 
correspond to the second sublattice (Fig. |l|). The choice 
of such a correspondence is not unique and it affects the 
form that the interaction takes after the transition from 
the two-leg ladder model to an equivalent 2D bipartite 
lattice. 




FIG. 1. 2D two-leg ladder: black and white circles com- 
prise the first and the second legs, respectively. Solid lines 
denote rungs. 



The kinetic term on the two-leg ladder includes next- 
neighbor hopping along the legs. For the equivalent 2D 
bipartite lattice, in general case, the hopping between re- 
mote sites will appear on the second sublattice. However, 
there existQ such a correspondence between the ladder 
and the bipartite lattice so that the kinetic term on the 
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bipartite lattice will include only next-neighbor hopping. 
Furthermore, the kinetic term will be SU (4)-invariant 
only if each of the sublattices is bipartite as well. Then 
the corresponding two-leg ladder will be bipartite too, 
and its rungs can be enumerated so that even and odd 
rungs will belong to different subladders. It is the bipar- 
tite 2D two-leg ladder that is considered in this paper. 

The studied model is of special interest not only be- 
cause of its probable relationship to high-T c cuprates, 
but also due to the possibility to .synthesize such ladders 
and explore them experimentally.!! 2 ! Therefore, one of the 
goals of the presented work is to derive the analytical 
expressions for the critical temperature of the transition, 
as well as for specific heat and for components of DC 
conductivity tensor in normal state. The calculation is 
performed in weak coupling limit and mean-field approx- 
imation, however, some of the results can be qualitatively 
extended to the strong coupling regime. 

In this regime there are six possible phases that corre- 
spond to 15 generators of the algebra and the Casimir 
operator. The values of the operators alternate their 
signs at the nearest neighbor rungs, although not all of 
the phases are density waves. In presence of symmetry 
breaking terms, only AF and dSC states remain. 

A peculiarity of the constructed model is that the in- 
teraction part of the Hamiltonian is being split into two 
terms with SO (4) x U (1) symmetry that renormalize in- 
dependently. Consequently, these terms correspond both 
to different conditions on phase transition and to the dif- 
ferent kind of the transition (AF or dSC). In particular, it 
will be shown that in mean-field approximation the pos- 
sibility of the dSC transition is determined by the value 
of chemical potential, but not by temperature, and there- 
fore, the system can be in the dSC state at arbitrary high 
temperatures. This property is the direct consequence of 
the presence of SU (4)-invariant interaction in the model 
that attracts the fermions of different "flavors" . 



II. MODEL 

In a two-leg ladder model each rung can have up to 
four fermions, differing in spin and "flavor" . This model 
is related to an equivalent 2D bipartite lattice if we 
assume that one flavor, c, corresponds to fermions on 
the first sublattice, while the other flavor, d, is a linear 
combination of the fermions on the second sublattice, 

(rj) = J2r' 4> ( r i - r ') c a (r'), where (r) = for x + y 
even and j is labeling the sites within the first sublattice. 
In the important caseS when cf>(r) = (4/-7r) (x 2 — y 2 ) 1 
for x + y odd, the operator Cko-dko- takes the form of 
the dSC order parameter Ak = <?kCko-C-ko-, where g^ = 
sign (cos k x — cosk y ) and the annihilation operators Cko- 
and c_k<r act on different sublattices. 

The definition above leads to a number of important 
consequences. First, we can regard <i-particles as well- 
defined fermions, since their creation and annihilation 



operators obey canonical commutation relations. Sec- 
ond, the total number of d-fermions is equal to the num- 
ber of fermions on the second sublattice. Third, in 2D 
case the contribution from the second sublattice to ki- 
netic energy takes the same form be it written in terms 
of initial c-operators or newly defined <i-operators. Thus, 
we conclude that d-fermions are the alternative represen- 
tation of the second sublattice of the 2D bipartite lattice. 

As it was mentioned in the introduction, we will as- 
sume that the two-leg ladder is bipartite. In order to 
simplify the construction of the Hamiltonian with explicit 
SU (4) symmetry, let us group the fermion operators on 
a rung into a 4-componcnt operator 



(cj T , c] i; (-l) j d jh (-l)Hi). (2.1) 



The terms that include only the scalar products of 
such operators are SU (4)-invariant. Those that involve 
the antisymmetric inner product ^j a E a ^^ reduce the 
symmetry of the group to Sp (4), or equivalently, SO (5). 
As we will see below, such terms bring about leg-to-leg 
hopping within a rung. 

Consider a model with the Hamiltonian 



He + Hchcm 



(2.2) 



Here the kinetic (hopping) and the scalar interaction 
terms are SU (4)-invariant: 
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where the "hypercharge" operator Yj = rij — n"- 1 ' — 



(2.4) 
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Ea4<r C J CT > and n 3 



(J) 



Thus, different flavors have opposite hypercharge. Also 
note that the kinetic term is invariant only globally, since 
it contains scalar products of the operators on different 
rungs, while the interaction term is locally invariant as 
well. 

The SU (4) — > SO (4) x U (1) breaking terms are chem- 
ical potential and Coulomb interaction. The latter can 
be regarded also as dSC - AF anisotropy, as it can be epu 
pressed in terms of the square of the local spin operator:!!!!! 
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(2.6) 



Here rij = n~- + nf\ One can also introduce the 
SU (4) — > SO (5) symmetry breaking terms so that their 
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combination with Eqs. ( pTB 
symmetry to SO (3) x U (1) 



2.6) will finally reduce the 
Such terms are the SO (5) 



"rotor" , or Casimir, operator and the inter-rung leg-to- 
leg hopping: 



Hmt - E Yl £2 j*b - E 

A j a<b j 



1 

2^ 



nf - 1 



nf - 1 



g(c) _ g(d) 



(2.7) 



^hop 



3<r 
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J J I J J 



h. c. 
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where £ ja b are the generators of SO (5) symmetry de- 
fined below in Eq. ( §Io| ), S (c) = ±cW, and = 



\Sad. In the absence of leg-to- leg hopping Eq. (2. 
the total hypercharge of the system ^ Yj is a conserved 
quantity. 

Finally, there exist a class of terms that break sym- 
metry to SU (3). They correspond to the presence of 
certain inhomogeneous magnetic field that takes differ- 
ent values on the sublattices. An example of such a term 
is -H< c ) • SH 

The free energy spectrum of both c and d-fermions is 
e (k) = —It (cos k x + cos k v ). The interaction term in 
Eq. (|2.2| ) can be written as a sum of 15 generators of 
SO (6) = SU (4) so that Hint becomes V ■ Hi n t,j with 



n 



mtj 



5 
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(2.9) 



Here M a {, are the generators of the matrix representation 
of SO (6) that acts on the space of 4-by-4 matrices by 
conjugation, C ab are the generators of the representation 
of SO (5) and N a is the corresponding superspin: 



£ ab =-¥M ab V fora,6=1...5, (2.10a) 
JV a = ^M 0a y for a = 1... 5. (2.10b) 

It is convenient to choose the following representation 
for M ab E 



M 0a 
M 04 
M ab 



<Ja 

ol 

— i<Ty 

ia„ 



1,2,3, 



cr„ 



--[M 0a ,M 0b ], a,b=1...5. 



The physical meaning of the components is M a = 
ReAq,ImAq) T , where antiferromagnetic 



order 



parameter m 



Aexp(— iQ ■ r), dSC order parameter A = iccr y d, Q 
(7r, 7r), and £ a (, incorporates rung spin operator S = 
S( c ) + S^, 7r-operator 71-t = —^aa y d^, and electric 
charge density Q — i (n( c ) + n^') — 1. In the absence of 
St/ (4) symmetry breaking, the components of C ab can 
also evolve into order parameters, such as ferromagnetic 
order parameter and 7r-wave superconducting order pa- 
rameter. Note that if we formally replace d- operators 
in Aq by c-operators, the result will coincide with the 
expression for the ry-operator that generates SO^L-useu- 
dospin symmetry in the standard Hubbard modelo 

All 15 generators of SU (4) algebra can be regarded as 
components of a superspin that transforms according to 
the adjoint representation of SU (4) . Although the su- 
perspin is not the order parameter, it is directly related 
to one. 

It is known, from the theory of the 2D standard Hub- 
bard modeltij that at half-filling below the transition 
point the density-wave state has lower ground energy 
than the spatially homogeneous state if the lattice is bi- 
partite. In such a state the components of the superspin 
will alternate the sign at even and odd rungs. Thus, if 
the ladder (or each of the sublattices of the equivalent bi- 
partite lattice) is bipartite, charge-density wave (CDW), 
spin-density wave (SDW), and dSC state are the only 
actually possible ordered phases at half-filling. In the 
pure SU (4) theoryli3 there are totally six phases that 
can be classified by the type of order and the parity with 
respect to the exchange of c and d-particles. The table 
below displays the components of the superspin that vary 
as cos (Q • r) for CDW and SDW states and the actual 
order parameters for the dSC states, according to such a 
classification: 



Order 


CDW 


SDW 


dSC 


odd parity 
even parity 


n {c) _ n (d) 
n (c) + n (d) _ 2 


gw _ gw 
g(c) + S (d) 


A, At 

7T, 7T 



(2.11) 



Note that the odd parity CDW phase takes place only 
when the coupling U is negative, while the rest only when 
it is positive. 

In the dSC states the components of the superspin vary 
as cos (Q • r), but due to the presence of (— l) 3 factor in 
Eq. (pTl]), these components (such as Aq) become nat- 
urally related to the quantities that are constant every- 



3 



where in the dSC state (such as A). The latter are the 
order parameters. In the case of CDW and SDW states, 
the order parameter is the amplitude of the variation, or 
the difference between the value of the superspin on even 
and odd rungs. 

There are also 16 eigenstates of Hint [Eq. (|2.9| )1 that 
can be labeled by the eigenvalues of M 2 — ^ a N 2 , 
£ 2 = EaEf,> £ L run S s P in component S z , charge 
density Q, and hypercharge Y (Table [jj). The ground 
state of Hint is 6-fold degenerate, consisting of the state 
|fi) = fcjdj — cjdj^ |0) and its five transformations 

by the components of the SO (5) superspin. Therefore, 
in the ground state of the total interaction term Hint, 
each rung is occupied only by c — d pairs. 

The degeneracy of the ground state of Hint at half- 
filling includes the contributions from the states, in which 
every rung of the two-leg ladder is occupied by one c — d 
pair, and the states with some of the rungs being fully 
occupied and some empty, or in other words, with four- 
particle c—d Cooper pairs. In the absence of c—d Cooper 
pairs, the degeneracy of the ground state is one of stan- 
dard Hubbard Hamiltonian, 2 , where N is a number 
of fermions on the lattice. Taking into account that k 
rungs can be empty and k rungs full, we can find the to- 
tal degeneracy and express it in terms of hypergeometric 
function F = iF\ (a, b; c; x): 



= 2 N F[ 



1 N N 



1 



xx-(T)(? 

k=0 v 7 v 



and for large values of N this expression can be roughly 
approximated as 5 N ^ 8 (8/N) 1 ^ 2 . 

In the strong coupling limit, the ground state is ap- 
proximately one of Hint an d we can derive the analog of 
the t — J model by computing the second-order correc- 
tion to the kinetic term (as zero and first orders vanish 
in the ground state). Using the identity M Q ^ • M 7< 5 = 
^asSfi-y — SapS-ys and taking into account that Yj ; = in 
the ground state of Hint, we find: 



H 



t-j 



(iJ) 



Mi ■ Mi 



tj -t, 



(2.12) 



where J = t 2 /U, Mj, ab = ^]M ab ^j, and t jtab = 
jMab^j. Note that in the given representation for 
M ab some of the components of tj vanish. The third 



term in Eq. (2.12) has a physical interpretation as pair 
hopping. 



III. CRITICAL TEMPERATURE 

In this section we will use the temperature Green func- 
tions technique to derive the expression for the crit- 
ical temperature of the phase tra nsitions in the sys- 
tem with the Hamiltonian Eq. (2.2). The bare fermion 



Green function is <S^1 (k, u) = 5 a p/ [icu n — Y a (e^ — fi)] , 
uj n = 2u (n + 1/2) /csT, where Greek indices run from 
1 to 4 and correspond to Cj, Cj, dt, and d\, re- 
spectively; Y a = +1 for a = 1,2, and —1 for 
a = 3,4. First, consider the case when g = 
0. The antisymmetrized bare interaction vertex cor- 
responds to the factor -T (7l ' 72;73 ' 74) (fci, A: 2 ; fc 3 , fc 4 ) = 

U (^7174^72 73 

has been used 



<5 7l73 <5 7274 ), where the notation k — (k, uj) 





s-q 



(a) (b) 

FIG. 2. (a) Self-energy in one-loop order, (b) Correction 
to interaction vertex. 

The chemical potential is fj, — /iq + Y a Y,( a > a ^ (0) , 
(a = 1, ...,4) where £( Q >' 3 ) (k, u n ) is the self-energy. The 
latter can be computed to one-loop order [Fig. 0(a)] and 
turns out to be a constant, T, (a '^ (k,w„) = (0). 
Therefore, in one-loop order the Green function does 
not change, since the correction to the single parti- 
cle energy and one to the chemical potential cancel, 
hence, the Fermi level remains Ef = fJ- — Y a ^ a ' a ^ (0) = 
fiQ. The corresponding density of states can be found 
from the expression for the free-particle energy = 
—2t (cos k x + cos k y ) and near the zero energy level £k = 
it diverges logarithmically: 

1 1 

8t 



p (e) ~ sign (e) 



(27T)" 



2t 



2.75 



lei < 2t. 



(3.1) 



Next we allow g to take a s mal l nonzero value. Then it 
is necessary to split H; n t Eq. (2.4) into a sum of two terms 
so that one of them will be similar to the Coulomb term 
Introduce two vectors PW and P^ 2 ' defined 

1 for a = 1,2 and for a = 3,4 



P (i) 



Eq. (f2_6|). 

as follows 

and Pi 2) = 1 - P y a 
can be represented as 



Then each term in Eq. 



where /; 



M 



M 



(i) 
j,ab 



(1) = A4 (1)2 
*tM«^, and M« 



1(1) 

l j,ab 



■(2) 



(2)2 h {2) - 2M (1) M {2) 
M„ h -PW. Thus, 



j,ab ' 



H 



1 



hit 



5 

3 a>b 

j a>b 



M {2) f 



V-V-„ . ,(1) M y 



<(2) 



(3.2) 
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up to an additive constant. In the diagrammatic cal- 
culations, the vertices 1^ and T^ 2 \ corresponding to 
the first and the second terms in Eq. (3.2) respectively, 
satisfy the identities 1^ = ToM^', z = 1,2, where 
{Ao B)^ 1,12 ' 13,1 ^ = J2 p A (7l ' 72;/3l '' ?2 )i3 ( ' 3l '' 32;73 ' 74) , 
and £/( z )(ti>T2;'»>T4) are antisymmetric with respect to the 
interchange 71 <-> 72 and 73 «-> 74 and are defined by the 
components: 

7 /(l)(12;12) _ 1 ? /(l)(34;34) _ I 

2' 2' 

7 /(2)(13;13) _ 1 7/(2)(14;14) _ 1 
2' 2' 



U 



(2)(23;23) 



I 7/(2)(24;24) _ \ 

2' " 2" 



The rest of the components of that remain unde- 
termined after antisymmetrization are zero. Then it 
follows that US Z ^ and 1^ have the following proper- 
ties: W« oZ/W = 0, T( 2 ) ol/W = o^) = 1<- z \ 

T( z ) o T( z ) = {U^/hfu^\ where f/M = U + g and 
(7< 2 ) = U. 

In order to calculate the critical temperature, we need 

(z) 

to evaluate the effect of the renormalization of 1 ab by 
the next-order corrections. There are three of them in 
the second order, and for the purpose of the calculation 
of interest, the primary contribution is given by [Fig. 
1(b)] 



J 



-1(71, 72:73,74) 

ab 



M) (k U k 2 ; S ~q,q) (q) ©£ 2 (s - q) . 



(3.3) 



r 



where s = ki + k2 = k3 + k4 and lo s = uj\ + u>2 = ^3+ W4 
are small. By considering simultaneously the corrections 
to the vertex (U/b)U^ and taking into acco unt that 
7Y int includes the sum over 15 generators [Eq. (2.9)], we 



Tz) 

find the expressions for the complete vertices, T b 



ab/ (i - y^U^K^ ,2 = 1,2, where 
d 2 q 



.(1) = ± 



rl „ tanh ^ 



,(2) 



2 J (2tt)' 
p(mo) • 



e (q) - Mo 



(3.4a) 
(3.4b) 



At critical temperature some of the complete vertices 
Xc^ diverge. It takes place when either of the following 
conditions is satisfied: 



(U + g)K<n = \/-, or 



The divergence of 2^ can be regarded as one of the 
Coulomb term Tic, even if initially the latter was ne- 
glected. It implies that the corresp ondin g transition is 
antiferromagnetic. Therefore, Eq. (3.5a) is the condi- 



(3.5a) 
(3.5b) 



tion on AF transition and, similarly, Eq. (3.5b) is one 



on d SC tr ansition. The integral in the right-hand side of 
Eq. ( 3.4a ) can be evaluated (see Appendix), which gives 
the value of the critical temperature of AF transition: 



/■■I," :~ 2/Dexp j -y jj— * 



2tD 



m 

2t 



3/8 



3\ 1/2 U + 9 



2t \ ksT, \M , 
U + g)' 2t 2t 



iMol k B T c 
2t 2t ' 



kBT c 
2t 



»1, 



(3.6a) 

(3.6b) 
(3.6c) 



where D = 72 1 / 4 //T 1 / 2 ps 0.387 and 7 w 0.577 is Eu- 
ler's constant. We see that if |/io| — 2tD 8 ^ 5 , the critical 
temperature can be as high as the difference between the 
Fermi level and the zero energy level £k = 0. For temper- 



atures high er th an 2t the critical temperature attains the 
value Eq. ( 3.6c ), although this assumes that U/2t 3> 1 
and therefore can not be regarded as a rigorous solution. 
However, the latter result allows one to make a quali- 
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tative conclusion that in the strong coupling limit the 
critical temp erature is of order of U. 

Eq. ( 3.5b| ) appears to be an equation on the critical 
value of chemical potential fiQ C so that the transition 
takes place at small values of /iq. This has to do with 
the fact that in one-loop order the density of states is di - 
verging at the zero energy level and this is why Eq. ( 3.5b| ) 
can be solved for /xo c for arbitrary weak interaction U. 
In next orders the density of states will likely to become 
finite and the peak will d ecrease with temperature, there- 
fore, the equation (3.5b) will evolve into one on the crit- 
ical temperature for the dSC transition. 

Finally, by-, explicitly diagonalizing the mean-field 
Hamiltonian,t2l one can find the spectrum of the exci- 
tations below critical temperature and the gap. For the 
AF transition, i.e., when one of the components of the su- 
perspin varies as (A/}) = — (5/8J7) N cos (Q • r), where 
Q = (7T,7r), j — 1, 2, or 3, the Hamiltonian eigenval- 

1/2 

ues are E± (k) = n ± (4 + AT 2 /4) ' (double de gener- 
ate). Similarly, for dSC transition (j = 4 or 5), there 



are four branches, E± (k) 



(/i ± e k ) 2 + A 2 /4 



and E { 1 ] (k) = -E. 



(i) 



1/2 



. (k). In the purely SU (4) casi 
/i = g = 0, the gap is equal N — (87r/e) t exp (— 5nt/U) 



IV. SPECIFIC HEAT AND ELECTRICAL 
CONDUCTIVITY 

The calculation in the previous section shows that in 
the mean-field approximation for a certain interval of the 
values of chemical potential the system is in dSC state at 
all temperatures. Certainly, if we included higher orders 
into our calculation, there would appear finite critical 
temperature. Although, strictly speaking, we have not 
shown existence of a transition to the normal phase as 
temperature increases, we can nonetheless consider the 
study of such a phase as a separate problem. 

The most interesting case is when the critical temper- 
ature for the dSC transition is of order or higher than 
hopping t. Then the normal state becomes automati- 
cally the high-temperature regime for the gas of fermions, 
which dramatically changes the temperature dependence 
of kinetic and thermodynamic quantities. 

The high-temperature limit of specific heat can be 
computed using the formulaaia 



1 



cv 



2 (k B TY 
I 

2 (k B T) 2 



d 2 k 



(27T) 

de e 2 p (e) , 



dw uj 2 Im 



G 



du> 



w=0 

(4.1) 



2h 2 a 



+ OC 



de K (e — (i) pi (s) 



(4.2) 



xA(w -Y a (e- n)) , 



Pi (e) 



d 2 k ( ctek 



(2tt) 2 \dk 



6(s- £k) , 



where the spectral weight A (u>) — — — Im (a; + i/2r)~ 1 , 
a is a lattice constant for the sublattices, and r is the 
relaxation time. 

In order to analyze the Hall effect, it is necessary to 
modify the kinetic term in the Hamiltonian by taking 
into account the presence of weak magnetic field: 

d 2 k 



a* - - 27r > V 



-te 



(27 

L a<T,k x ,fc B -() Q L ao'k 

— It cos kyC aa j l c a(7 ] i j , 



(4.3) 



where b a — Y a b, b — \e\a 2 B z /hc, and we have assumed 
that the carriers are negative-charged. The xy compo- 
nent of the conductivity tensor is 



'xy 



P2 (e) 



2 \e\aB z 
3^ 4 c 

d 2 k ( de 



de 2tK (e - ft) p 2 (e) , (4.4) 



8 2 e 



{2^f\dkJ dk 2 /^ 



and the Hall coefficient Rh — cr xy /a 2 x . In the high- 
temperature limit kgT 3> t, the derivative of the Fermi 
function dfp/duj ~ l/4fcsT and consequently, in the 
same approximation as in the previous section, we find 



h 2 a ksT' 



3^2 



'xy 



Rh 



8 \e\ aB z t A r 

3H 4 c ksT' 
8a 3 B z k B T 
3 lei c t 



(4.5a) 
(4.5b) 
(4.5c) 



As we can see, the Hall coefficient increases linearly with 
temperature. On the contrary, in the low-temperature 
limit, we recover the usual linear dependence of specific 
heat and weak parabolic dependence of DC resistivity on 
temperature as well as no dependence of the Hall coeffi- 
cient on temperature. 



V. CONCLUSION 



where Gr {uj, k) is the retarded Green function. Thus, 
c v oc {t/TfioT k B T > t. 

The DC conductivity can be derived from Kubo for- 
mula. In the large relaxation time approximation rt 3> 1, 



I have studied the model on a 2D bipartite lattice, 
which is equivalent to the two-leg ladder, that includes in- 
teraction of fermions from opposite sublattices. The cor- 
responding term in the Hamiltonian is SU (4)-invariant. 



G 



The symmetry breaking factors include chemical poten- 
tial and Coulomb interaction. Fermions that belong 
to different sublattices have opposite "hypercharge" , a 
symmetry-related quantum number. In the absence of 
leg-to-leg hopping the total hypercharge of the system 
conserves. 

The physical meaning of the components of the order 
parameter is slightly different from those given in Ref. @, 
primarily in that the ground state of the "antiferromag- 
netic" phase is actually a density wave with varying rung 
Neel vector. In the strong SU (A) interaction limit, the 
degeneracy of the ground state is higher than one of the 
standard Hubbard model and the Hamiltonian resembles 
one oft — J model, but also includes a pair hoppin g te rm 
with operators similar to those introduced in Ref. |l7]. 

The condition on the phase transition has been evalu- 
ated in one- loop order in weak interaction limit. It leads 
to the equation on critical temperature for the AF transi- 
tion. However, for the dSC transition, it depends on the 
density of states at the Fermi level, but not on tempera- 
ture. In the latter case, the dependence on temperature 
will probably appear in the next orders of perturbation, 
but the resulting value of critical temperature will be 
still abnormally high. Furthermore, if we formally con- 
sider the strong coupling limit in the derived formulae, 
the condition on dSC transition will be always satisfied, 
except when the Fermi surface lies in the area where the 
density of states is smaller than a threshold value of order 
of inverse coupling coefficient. Since the density of states 
will have a drop at the zero energy level, the AF tran- 
sition will occur near to half-filling and dSC transition 
away from half-filling. When chemical potential Mo = 0, 
no transition at finite temperatures takes place due to 
the Mermin - Wagner theorem. 

The fact that the critical temperature of the dSC tran- 
sition can be of order or higher than hopping implies that 
in the normal state the gas of fermions will no longer 
resemble the usual zero-temperature Fermi liquid with 
discontinuity of the distribution function at the Fermi 
surface. The distribution function will almost linearly de- 



crease with energy, which will result in abnormal temper- 
ature dependence of all experimentally measured charac- 
teristics. In particular, specific heat will be inverse pro- 
portional to the square of the temperature, and the DC 
electrical resistivity and the Hall coefficient will be lin- 
early proportional to the temperature. Such dependence 
agrees with experimental data on high-T c cuprates.ES 
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APPENDIX: 

In this Appendix the integral that appears in the right- 
hand side of Eq. (3.4a) is evaluated: 



2tt 
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tanh 



f e(k)-Mo "\ 
V 2fe B T c J 
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(Al) 



First, we make a substitution fc+ = (k x + k y ) /2, k— = 
(k x — k y ) /2 and expand the energy in terms of k + about 
the Fermi level e (k) = /j,q: 



e (k) - fio 

x k+ — arccos 



— y (At cos k- ) — Mo 

Mo 
At cos k- 



(A2) 



In the limit of Mo/2i — > 0, the integral I diverges loga- 
rithmically. The internal integral can be taken by parts, 
which results in a logarithmic part and a convergent in- 
tegral. In the latter the limits can be replaced by ±oo. 
Then there will be a region of integration at k x = 0, 
k y = 7r/2 that will not be covered and a symmetric region 
that will be covered twice. However, the expression in the 
integral takes the same value in both regions, therefore 
the result remains unchanged and after the limit replace- 
ment no corrections will be necessary: 
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dk^ 
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7 I 



(A3) 



Here Euler's constant 7 » 0.577 and mo is assumed to 
be positive. Furthermore, Eq. (Al) does not depend on 
the sign of Mo, thus, we can replace Mo by I Mo I- The 
asymptotic expansion at |mo| /2i — ► is 



AtxH 
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TABLE I. Classification of the eigenstates of Timt- 
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